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Abstract 

This  paper  presents  the  work  on  a  simple,  natural  gas-fed,  hybrid  solid  oxide  fuel  cell-gas  turbine  (SOFC-GT)  power-generation  system. 
The  system  consists  of  an  internal-reforming  SOFC  (IRSOFC)  stack,  a  combustor,  a  GT,  two  compressors  and  three  recuperators.  Two  case 
studies  are  conducted  with  particular  attention  on  the  effects  of  operating  pressure  and  fuel  flow-rate  on  the  performance  of  the  components 
and  overall  system.  Results  show  that  an  internal-reforming  hybrid  SOFC-GT  system  can  achieve  an  electrical  efficiency  of  more  than  60% 
and  a  system  efficiency  (including  waste  heat  recovery  for  co-generation)  of  more  than  80%.  It  is  also  found  that  increasing  the  operating 
pressure  will  improve  the  system  efficiency,  whereas  increasing  the  fuel  flow-rate  (while  keeping  the  fuel  utilisation  rate  unchanged)  causes 
the  system  efficiency  to  decrease.  In  the  latter  case,  the  increase  in  system  fuel  consumption  is  relatively  higher  which  removes  the  benefit  of 
increase  in  SOFC  stack  and  turbine  power  output.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

The  demand  for  power-generation  systems  of  high  effi¬ 
ciency  and  with  low  emission  is  of  increasing  importance. 
Fuel  cells,  as  alternatives  to  conventional  energy-conversion 
systems,  have  the  prospect  for  exploiting  fossil  fuels  more 
benignly  and  efficiently.  The  various  types  of  fuel  cells 
symbolise  quite  different  technologies,  with  no  clear  winner 
to-date  [1].  Nevertheless,  the  high-temperature  solid  oxide 
fuel  cell  ( SOFC)  appears  better  suited  to  power  generation  in 
a  hydrocarbon  fuel  economy.  Moreover,  due  to  the  syner¬ 
gistic  effects,  hybrid  SOFC  and  gas  turbine  (GT)  technology 
has  been  identified  to  be  a  superior  power-generation  tech¬ 
nology  to  many  other  options.  Predicted  results  shown  that 
the  integration  into  a  more  complex  GT  cycle  of  a  pres¬ 
surised  SOFC  which  is  multi-stage,  inter-cooled,  reheated 
and  recuperated,  can  result  in  an  electrical  efficiency  of 
70%  or  higher  (net  ac  per  lower  heating  value  (LHV )  of  fuel) 
can  be  expected  [2—4],  For  small-scale,  electrical  power- 
generation  systems  consisting  of  advanced  microturbines 
and  SOFC  stacks,  with  a  power  rating  of  250  kW,  an  elec¬ 
trical  efficiency  of  65%  (net  ac  per  LHV  of  fuel)  can  be 
expected  [5],  Other  advanced  system  studies,  which  include 
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thermo-economic  assessment  of  internal-reforming  SOFC 
(IRSOFC)-GT  combined  cycles  and  take  into  account  seve¬ 
ral  technological  constraints,  have  also  been  reported  in 
detail  [6], 

In  this  investigation,  attention  is  focused  on  the  develop¬ 
ment  of  a  model  for  a  simple  natural  gas-fed,  hybrid  SOFC- 
GT  power-generation  system.  The  model  is  based  on  the  first 
law  of  thermodynamics.  The  IRSOFC  stack,  which  is  the 
heart  of  the  system,  interacts  with  other  system  components 
and  causes  a  complicated  interdependency  between  them. 
Thus,  it  is  useful  to  consider  all  performance-related  para¬ 
meters  as  a  function  of  the  thermodynamic  states,  e.g.  to 
relate  the  exchange  current  density  of  the  fuel  cell  and 
resistivities  of  cell  components  to  the  temperature  of  the 
SOFC  stack.  A  general-purpose  mathematical  software  tool, 
MATLAB,  is  adopted  to  develop  the  simulation  code. 

2.  System  configuration  and  description 

The  configuration  of  the  hybrid  SOFC-GT  power-gen¬ 
eration  system  is  shown  in  Fig.  1.  The  system  consists  of 
an  IRSOFC  stack,  a  combustor,  a  GT  and  a  power  turbine 
(PT),  two  compressors,  and  three  recuperators.  An  inverter 
(not  shown)  with  a  typical  efficiency  of  96.5%  is  considered. 
The  mechanical  work  produced  by  the  GT  is  used  to  drive 
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Nomenclature 

A  constant 

B  constant 

cp  specific  heat  at  constant  pressure  (J/mol  K) 

E  cell  potential  (V) 

£°  reversible  cell  potential  (V) 

F  Faraday  constant  (96  487  C/mol) 

G  Gibbs  free  energy  (J/mol) 

h  enthalpy  (kJ/mol) 

i  current  density  (mA  cm-2) 

;'o  exchange  current  density  (mA  cm~2) 

iL  limiting  current  density  (mA  cm"  2) 

Kp  equilibrium  constant 

iij  molar  flow-rate  of  specie  i  (kmolh-1) 

p  pressure  (bar) 

Q  heat  transfer  (kJ) 

Qrxn  heat  due  to  reactions  (kJ) 

Qsurr  heat  loss  to  surroundings  (kJ) 

R  universal  gas  constant  (8.314  J  mol  K) 

R,  ohmic  resistance  (Q) 

S  entropy  (J  mol  K) 

T  temperature  (K) 

Uf  fuel  utilisation 

Z  number  of  electrons  participating  in  the 

reactions 

Greek  letters 

P  transfer  coefficient 

8  thickness  (cm) 

t] Act  activation  polarisation  (V) 

t/conc  concentration  polarisation  (V) 

t/ohm  ohmic  overpotential  (V) 

p  resistivity  (Q  m) 

Subscripts 
a  anode 

c  cathode 

r  reforming  reaction 

re  reversible 

s  shifting  reaction 


the  two  compressors,  and  the  sole  mechanical  power  available 
for  electrical  power  generation  is  derived  from  the  PT. 

Compressed  fuel  and  air  are  preheated  at  their  respective 
recuperators  prior  to  entering  the  IRSOFC  stack.  Natural 
gas  and  oxygen  are  channelled  through  the  anode  and 
cathode  compartments,  respectively.  Natural  gas  is  intern¬ 
ally  reformed  in  the  anode  compartment  and  hydrogen-rich 
products  are  produced.  The  electrochemical  reaction  occurs 
at  the  three-phase  boundaries  (TPB)  of  both  electrodes,  and 
produces  ionic  flow  through  the  electrolyte  and  electron  flow 
across  the  electrodes.  Electrical  work  (energy)  is  hence 
produced  together  with  heat  generation.  The  generated  heat 
is  partly  dissipated  to  the  environment,  partly  used  to  reform 


the  natural  gas,  and  partly  used  to  heat  up  the  feedstock  and 
effluent  gases.  The  high-temperature  effluent  gases  of  the 
IRSOFC  stack,  i.e.  unutilised  reformed  natural  gas  and 
depleted  air,  are  channelled  to  a  combustor,  where  residual 
fuel  (hydrogen,  methane,  carbon  monoxide)  reacts  with 
the  excess  air.  The  amount  of  residual  fuel  in  the  exit  of 
the  SOFC  stack  is  determined  by  the  operating  conditions 
of  the  IRSOFC. 

The  combusted  fuel-air  mixture  then  flows  through  the 
GT  whereby  mechanical  power  is  produced.  The  produced 
power  is  consumed  by  the  compressors  which  being  the  fuel 
and  air  to  their  required  pressures.  Subsequent  expansion 
through  the  PT  produces  additional  mechanical  power  which 
is  used  to  generate  the  electrical  power.  The  effluent  gases 
exiting  from  the  PT  pass  through  a  series  of  three  recup¬ 
erators  where  preheating  of  fuel  and  air  occurs.  Recuperator 
3  is  used  for  general  heating  purpose  similar  to  that  in  a  co¬ 
generation  plant.  For  example,  it  heats  liquid  water  from  25 
to  90  °C  for  industrial  processes. 


3.  System  modelling 

3.1.  Internal  reformer  model 

For  a  natural  gas  fuelled,  hybrid  SOFC-GT  system,  either 
internal  or  external  reforming  is  needed.  To  reduce  the 
expense  of  an  external  reformer  and  to  provide  additional 
‘cooling’  for  the  SOFC  stack,  the  use  of  an  internal  reformer 
is  the  better  choice. 

In  the  reforming  process,  steam  is  needed  and  can  be 
obtained  by  vapourising  the  feed-water  using  the  waste  heat 
from  the  SOFC  stack.  The  reaction  mechanisms  for  the 
internal  reforming  of  SOFC  stack  are: 


x  =  CH4  +  H20 

->  CO  +  3H2  (reforming) 

(1) 

y  =  CO  +  H20  - 

>  C02  +  H2  (shifting) 

(2) 

z  =  H2  +  i(02)- 

->  H20  (electrochemical) 

(3) 

where  x,  y  and  z  indicate  the  respective  molar  flow-rate  in  the 
reaction.  The  equilibrium  constants  of  reforming  and  shift¬ 
ing  processes  are  temperature  dependent  and  can  be 
obtained  from  the  following  equation: 

log  Kp  =  AT 4  +  BT3  +  CT2  +DT  +  E  (4) 

where  constants  A,  B,  C,  D  and  E  are  listed  in  Table  1  [6,7]. 


Table  1 


Values  of  equilibrium  constants  of  reforming  and  shifting  processes 


Reforming 

Shifting 

A 

-2.63121  x 

10~n 

5.47301  x 

10~12 

B 

1.24065  x 

io-7 

-2.57479  x 

10~8 

C 

-2.25232  x 

10~4 

4.63742  x 

10~5 

D 

1.95028  x 

10^' 

-3.91500  x 

10~2 

E 

-6.61395  x 

101 

1.32097  x 

10* 
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Fig.  1.  Schematic  of  recuperated  internal-reforming  hybrid  SOFC-GT  system. 


Assuming  that  the  reforming  and  shifting  reactions  are 
always  in  equilibrium,  the  equilibrium  constants  can  also  be 
calculated  from  the  partial  pressures  of  the  reactants  and 
products,  i.e. 


„  _  pLpco 

Apr 

Pch4Ph2o 

(reforming) 

^  _  Ph2Pco2 

^ps  — 

PcoPh2o 

(shifting) 

Assuming  that  x,  y  and  z  are  the  molar  flow-rates  of  CH4, 
CO  and  H2,  respectively,  participating  in  the  reactions,  then: 


Qs  =  y{hco2  +  h  h2  —  he  o  —  /?h2o)  (U) 

Assuming  that  the  heat  generated  from  the  electrochemi¬ 
cal  process  of  the  fuel  cell  stack  is  (9ixn,  the  total  heat 
available  for  transfer  from  the  stack  will  be: 

Q  =  Q™-Qr-Qs  (12) 

3.2.  SOFC  model 

The  fuel  cell  model  developed  in  this  study  is  based  on  a 
tubular  design  and  is  fed  by  the  natural  gas,  though  it  can  be 


([CO]1  +  x  -  y)/p»;ot  +  2x)([H2]’  +  3x  +  y-  z)/(n\ot  +  2x)  2 
([CH4]‘  - x)/ (n[ot  +  2x)([H20]‘  -x-y  +  z)/(n\ot  +  2x)  “U 


([CQ2]'  +  y)/(»;ot  +  2x)([H2]‘  +  3x  +  y  -  z)/ (»[ot  +  2x) 
([CO]1  +  x - y)/(n\ot  +  2x)([H20]‘  -x-y  +  z)/ «t  +  2x) 


z=Uf(3x  +  y)  (9) 

where  superscript  i  is  inlet;  U[  fuel  utilisation  rate;  «[ot  total 
inlet  molar  flow-rate  (or  iiNq  +  /ih,o)- 

When  the  temperature  is  known,  the  equilibrium  con¬ 
stants  can  be  calculated  from  Eq.  (4)  and  the  unknown  x,  y 
and  z  are  determined  by  solving  the  simultaneous  Eqs.  (7)- 
(9)  at  the  given  fuel  utilisation  rate  and  inlet  conditions  of  the 
flow. 

Both  the  reforming  and  shifting  reactions  are  endother¬ 
mic;  the  heat  needed  for  each  reaction  can  be  calculated 
from  Eqs.  (10)  and  (11),  respectively,  i.e. 

Qr  =  x(hco  +  3hH2  —  hHlo  —  hcH4 )  (10) 


customised  to  accept  different  types  of  fuel  cell  designs  and 
feedstock.  The  geometric  and  performance  related  data  of 
the  SOFC  are  based  on  the  design  developed  by  Siemens- 
Westinghouse  [8].  To  simplify  the  study,  this  model  focuses 
on  the  thermodynamic  aspect  and  the  associated  electroche¬ 
mical  processes  of  the  cell  operation.  Treating  all  the  cells  as 
identical,  the  performance  of  the  stack  can  be  estimated. 
Hence,  knowing  the  amount  of  fuel  consumed  by  the  stack, 
the  total  current  produced  by  the  stack  can  be  calculated; 
while  the  cell  voltage  can  be  determined  based  on  informa¬ 
tion  for  the  stack  temperature  and  the  fuel  utilisation  rate. 
With  known  cell  voltage  and  the  total  current,  the  electrical 
power  of  the  stack  is  calculated. 
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For  a  H2-O2  reaction  in  a  SOFC  with  known  partial 
pressures  of  reactants  and  products,  the  ideal  reversible  cell 
potential  can  be  calculated  from  the  Nernst  equation,  i.e. 


£°  +  AE 


-A  G°  -AG 
nF  nF 


-A  G°  RT  Pn2Po2 
IF  IF  ph2o 
(13) 


There  are  several  sources,  however,  which  contribute  to  the 
irreversible  losses  in  an  operating  fuel  cell.  These  losses, 
usually  termed  polarisation,  overpotential  or  overvoltage  (;/), 
originate  primarily  from  three  sources,  namely:  (i)  activation 
overpotential  0/Act);  (ii)  ohmic  overpotential  (i] ohm);  and  (iii) 
concentration  overpotential  ( f/conc )•  Hence,  the  actual  poten¬ 
tial  of  an  operating  fuel  cell  will  be  in  the  form  of: 

^  The  H/ Act  T  hohm  3”  ^Conc)  (14) 

where 


V  Act  =  f  1 


or 


io(T) 


(1 


'lOhm  ~ 


i 


AjeP‘,TSj 

Asj 


^1  Cone 


where  />  is  the  transfer  coefficient;  i0  is  the  apparent 
exchange  current  density;  Rj  is  the  total  cell  resistance, 
which  includes  electronic  and  ionic  resistance;  Aj  and  P>t  are 
constants  specific  to  the  materials  used;  8j  is  the  thickness  of 
the  respective  cell  component;  zL  is  the  limiting  current 
density. 


3.3.  SOFC  stack  model 


The  SOFC  stack  consists  of  a  series  of  many  single  cells  in 
order  to  increase  the  voltage  output  and  hence  the  electrical 
power  for  practical  applications.  Since  the  SOFC  stack  inter¬ 
faces  with  many  other  system  components,  depending  on 
what  the  system  configuration  is,  its  performance  will  be 
closely  related  to  the  thermodynamic  states  of  the  components 
to  which  it  is  connected. 

To  calculate  the  stack  temperature,  three  heat  sources/ 
sinks  should  be  determined;  namely:  (i)  heat  generated  due 
to  the  electrochemical  reaction  in  the  SOFC  stack;  (ii)  heat 
consumed  by  internal  reforming;  (iii)  heat  loss  to  the 
surroundings. 

The  heat  generated  by  electrochemical  reaction  comes 
from  two  main  sources.  One  is  due  to  the  reversible  reaction, 
while  the  other  is  due  to  the  process  irreversibility  associated 
with  cell  polarisation.  The  reaction  heat  takes  the  form  of: 

Qaa  =  +  T  As  (!5) 


AS-(s°„,o  jS°o,)  2“('pCj 

where  tj  and  AS  are  both  temperature  dependent. 

Under  steady-state  operation  of  a  SOFC  stack,  the  heat 
generated  by  the  electrochemical  reaction  will  provide  the 
heat  for  internal  reforming  and  will  heat  up  the  mixture  of 
products  and  non-reacted  residual  gases.  By  applying  the 
steady  flow  energy  equation  and  assuming  that  the  amount  of 
heat  dissipated  to  the  surroundings  is  Qsurt  and  the  effluent 
gases  from  the  anode  and  cathode  exits  are  well  mixed  and 
arrive  at  temperature  T2,  then: 

Q  =  Gsurr  +  A/zci  +  A/?c2  +  Aha\  +  A/za2  (16) 

where  A/zci ,  A/;ai  are  the  enthalpy  changes  of  reactants  at  the 
cathode  and  anode,  respectively;  Ahc2,  Aha2  are  the  enthalpy 
changes  of  products  at  the  cathode  and  anode,  respectively. 
Given 


fT2 

Ah  =  /  cp  dr 


(17) 


and  combining  Eqs.  (12),  (16)  and  (17),  then: 


G'  =  Grxn  -  Gr  -  Gs  -  Gsurr 

rT 2  \  /  r  'i'2 

=  1  «i  /  cpdT  +  n2  /  cp  dT 


I A  a 


«3 


r*  72 


CpdT 


rT2 


ti4  /  Cp  dT 
3  Tic 


(18) 


where  n\,  n2,  n3  and  n4  are  the  corresponding  molar  flow-rates 
of  the  gases. 

Since  specific  heat  cp  is  also  a  function  of  temperature,  it  is 
fitted  to  a  second  order  polynomial  of  the  temperature  func¬ 
tion  and  is  integrated  across  a  temperature  span  of  7)  and  T2. 
In  this  study,  the  SOFC  stack  is  assumed  to  be  in  thermal 
equilibrium  with  the  effluent  gases.  Thus,  the  SOFC  stack 
effluent  gas  temperature  T2  can  be  determined  by  the  iterative 
computation.  The  iterations  continue  until  the  error  between 
two  successive  calculations  is  less  than  the  assigned  tolerance. 
The  error  to  meet  the  convergence  criteria  is  given  as: 


Gerror 


6'  -  (Grxn  Gr-  Gs  ~  6suff) 
Gixn  —  Gi  —  Gs  -  G  SUIT 


<  1% 


(19) 


3.4.  Other  system  component  models 

Other  system  components  that  also  play  crucial  roles  in 
determining  the  overall  performance  of  the  plant  include  the 
compressors,  turbines,  recuperators  and  combustor.  In  this 
study,  it  is  assumed  that  all  system  components  are  working 
at  their  respective  designed  conditions  under  steady-state 
operation.  A  set  of  operating  parameters  and  the  assumed 
efficiencies/effectiveness  of  these  system  components  are 
given  in  Table  2.  In  addition,  the  initial  states  of  the  working 
fluids  have  to  be  assumed,  in  particular  for  the  recuperators 
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Table  2 


Setting/assumed  values  of  performance-related  parameters 


Parameter 

Setting/assumed  value 

Compressor  isothermal  efficiency 

85% 

Recuperator  effectiveness 

0.80 

Recuperator  pressure  drop 

1% 

GT  mechanical  efficiency 

99% 

PT  isentropic  efficiency 

85% 

Inverter  efficiency 

96.5% 

Generator  efficiency 

97% 

SOFC  length 

50  cm 

SOFC  surface  area 

270  cm2 

Combustor  isentropic  efficiency 

95% 

Limiting  current  density 

350  mA  cm~2 

Electrolyte  resistivity  constants 

A  =  0.00294,  B  =  10  350 

Electrolyte  thickness 

0.004  cm 

Interconnect  resistivity  constants 

A  =  0.1256,  B  =  4690 

Interconnect  thickness 

0.01  cm 

Cathode  resistivity  constants 

A  =  0.00811,  B  =  600 

Cathode  thickness 

0.200  cm 

Anode  resistivity  constants 

A  =  0.00298,  B  =  -1392 

Anode  thickness 

0.015  cm 

Stack  pressure  drop 

2% 

Area  ratio  of  cathode  to  anode 

10 

as  the  hot  fluid  temperatures  (for  example,  the  flow  tem¬ 
perature  at  the  exit  of  the  PT)  are  not  available  at  the 
beginning  of  the  simulation.  These  assumed  initial  values, 
however,  will  converge  after  a  few  loops  of  iteration  until  the 
converging  criteria  are  met.  Computation  of  the  heat- 
exchange  between  the  hot  and  cold  fluids  in  the  recuperators, 
which  are  counter-flow  heat  exchangers,  is  based  on  the 
effectiveness-number  of  transfer  units  (e-NTU)  method. 

4.  Overpotentials  in  SOFC 

4.1.  Activation  overpotential 

The  activation  overpotential  at  different  operating  tem¬ 
peratures  using  the  Butler- Volmer  equation  is  shown  in 
Fig.  2.  The  results  show  that  the  activation  overpotential 
increases  with  increase  in  current  density,  but  decreases  with 
increase  in  temperature.  It  is  obvious  that  the  activation 
overpotential  is  zero  at  zero  current  density,  which  can  never 
be  the  case  if  the  Tafel  equation  is  used.  The  latter  yields  an 
unreasonable  activation  overpotential  of  —  ((RT)/ (4 fSF) ) In  io 
at  zero  current  density. 

4.2.  Ohmic  overpotential 

The  ohmic  resistance  of  the  various  cell  components  at 
different  temperatures  is  shown  in  Fig.  3.  It  is  obvious  that, 
based  on  data  for  the  Siemens-Westinghouse  SOFC,  the 
interconnector  contributes  the  most  to  the  overall  ohmic 
resistance,  and  is  followed  by  the  electrolyte.  The  ohmic 
resistance  of  the  electrodes,  by  comparison,  contributes 
negligible  overall  ohmic  resistance  in  the  SOFC.  Note  that. 


the  resistance  of  each  component  is  determined  not  only  by 
the  respective  resistivity  but  also  by  the  thickness.  It  can  be 
deduced  that  improved  performance  of  the  SOFC  can  be 
expected  if  a  cell  component  of  low  resistivity  is  used  as 
the  support/substrate  of  SOFC.  The  resistance  increases 
exponentially  with  decrease  in  temperature  and  this  result 
in  deteriorated  fuel  cell  performance  at  undesirable  low- 
temperature  operation. 

4.3.  Concentration  overpotential 

The  effect  of  the  temperature  on  the  concentration  over- 
potential  is  given  in  Fig.  4.  It  can  be  seen  that  there  is  a 
significant  increase  in  concentration  overpotential  when  the 
operation  of  SOFC  approaches  the  limiting  current  density 
of  345  mA  cm-2.  More  work  can  be  done  in  the  future  to 


Electrolyte,  interconnect,  cathode  and  anode  resistance 


Fig.  3.  Effect  of  temperature  on  ohmic  resistance. 
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Fig.  4.  Effect  of  temperature  on  concentration  overpotential. 

relate  the  effect  of  the  fluid  properties  on  the  limiting  current 
density  or  to  use  a  more  generalised  polarisation  model  [9]. 
Normally,  fuel  cells  are  designed  to  operate  slightly  lower 
than  the  current  density  that  corresponds  to  the  peak  power 
for  a  good  compromise  between  the  cell  efficiency  and 
operational  stability.  Hence,  under  normal  operating  condi¬ 
tions  the  present  model  is  still  applicable,  as  the  concentra¬ 
tion  overpotential  is  not  the  rate-limiting  factor. 

5.  Hybrid  IRSOFC-GT  system  simulation 

5.1.  System  performance  simulation  results 

The  pressure,  temperature,  gas  flow-rate  and  composition 
at  each  node  of  the  system  (refer  to  Fig.  1)  are  listed  in 
Tables  3  and  4.  In  this  particular  case  study,  the  pressure  ratio 
of  both  compressors  is  set  to  9,  the  air  flow-rate  to 


Table  4 

Gas  composition  at  different  nodes  of  the  system 


Gas  composition 


h2  o2 

n2 

H20 

n 

o 

ch4 

CO 

Stack  outlet  (anode)  (%) 
15.69  0 

0.32 

62.75 

16.19 

0.06 

4.99 

Stack  outlet  (cathode)  (%) 

0  9.36 

90.64 

0 

0 

0 

0 

Combustor  outlet 

0  3 

64 

26 

7 

0 

0 

170  kmol  h_1,  the  fuel  flow-rate  to  15  kmol  h_1,  and  the 
fuel  utilisation  rate  to  0.8.  The  total  number  of  tubular 
SOFCs  used  in  all  cases  is  60  000. 

Some  results  relating  to  cell  and  system  performance,  as 
well  as  the  conditions  at  the  exhaust  (node  13),  are  extracted 
from  the  simulation  and  listed  in  Table  5.  The  cells  operate  at 
power  density  of  104.5  mW  cm”2  at  current  density  of 
141.6  mA  cm-2,  which  is  far  away  from  the  limiting  current 
density  of  345  mA  cm  2.  The  waste  heat  recovered  from 
recuperator  3  is  731  kW  with  an  exit  flow  temperature  and  a 
pressure  of  393  K  and  1.34  bar,  respectively.  Note  that,  the 
net  electrical  power  output  (ac)  of  the  plant  is  2105.25  kW, 
of  which  381  kW  comes  from  the  PT.  The  latter  contributes 
about  18%  of  the  total  electrical  power  output.  The  electrical 
and  total  efficiencies  of  the  plant  are  62.2  and  83.8%, 
respectively.  The  total  (system)  efficiency  is  defined  as 
the  ratio  of  the  useful  energy  output  to  the  LHV  of  the  fuel; 
the  former  includes  the  electrical  energy  produced  by  the 
SOFC  stack  and  the  PT  (in  terms  of  ac  power  output),  and 
the  heat  recovery  at  recuperator  3.  Excluding  the  heat 
recovery  in  the  calculation  of  useful  energy  output,  the 
efficiency  so  obtained  is  defined  as  the  electrical  efficiency. 

5.2.  System  performance  under  different  operating 
pressures 


Table  3 

Fluid  properties  for  recuperated  hybrid  IRSOFC-GT  power  system 


Pressure 

Temperature 

Gas 

Gas  composition  (%) 

(bar) 

(K) 

flow-rate 
(kmol  h-1) 

1 

1 

300 

170 

N2:  79,  02:  21 

2 

9 

562 

168.3 

N2:  79,  02:  21 

3 

8.3 

960 

166.6 

N2:  79,  02:  21 

4 

1 

300 

15 

CH4:  97,  C02:  1.5,  N2:  1.5 

5 

9 

562 

14.85 

CH4:  97,  C02:  1.5,  N2:  1.5 

6 

8.55 

746 

14.7 

CH4:  97,  C02:  1.5,  N2:  1.5 

7 

7.64 

1166 

145.2 

Refer  to  Table  4,  cathode 

8 

7.65 

1166 

68.18 

Refer  to  Table  4,  anode 

9 

6.46 

1466 

206.34 

Refer  to  Table  4,  combustor 

10 

1.41 

1061 

204.3 

Refer  to  Table  4,  combustor 

11 

1.39 

791.5 

202.3 

Refer  to  Table  4,  combustor 

12 

1.36 

773.4 

200.2 

Refer  to  Table  4,  combustor 

13 

1.34 

393 

198.2 

Refer  to  Table  4,  combustor 

The  effect  of  operating  pressure  on  stack  temperature,  cell 
voltage,  stack  and  PT  ac  power  outputs,  waste  heat  recovery 

Table  5 

Cell  and  system  performance  related  results  from  simulation 


Cell  voltage 

0.738  V 

Cell  current  density 

141.6  mA  cm 

Gas  temperature  at  node  13 

393  K 

Gas  pressure  at  node  13 

1.34  bar 

Stack  temperature 

1166  K 

Compressor  work  (air  side) 

301  kW 

Compressor  work  (fuel  side) 

26.6  kW 

Net  PT  ac  power  output 

381  kW 

GT  inlet  temperature 

1466  K 

GT  inlet  pressure 

6.46  bar 

Waste  heat  recovery 

731  kW 

Electrical  efficiency  (LHV) 

62.2% 

Total  efficiency  (LHV) 

83.8% 
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Fig.  5.  Effect  of  flow  pressure  on  stack  temperature. 


and  system  efficiency  is  shown  in  Figs.  5-10,  respectively.  In 
this  case  study,  the  fuel  and  the  air  flow-rates  are  fixed  at  15 
and  170  kmol  h_1,  respectively,  with  a  80%  fuel  utilisation 
rate.  For  simplicity,  the  pressures  of  the  fuel  and  air  are  set  to 
be  the  same,  though  they  can  be  different  under  typical 
operational  demand.  Note  that,  the  increase  in  fuel  and  air 
pressures  through  their  respective  compressors  is  associated 
with  the  increase  in  flow  temperatures.  This  subsequently 
affects  the  operating  temperature  of  the  fuel  cell  through  the 
thermal  effect,  and  hence  the  overpotentials  and  heat  gen¬ 
eration.  This  means  that  the  higher  is  the  feedstock  tem¬ 
perature,  the  higher  will  be  the  stack  temperature.  A  higher 
stack  temperature  means,  however,  that  the  fuel  cell  operates 
more  efficiently  and  less  waste  heat  will  be  generated.  The 


Fig.  7.  Effect  of  flow  pressure  on  stack  ac  power  output. 


combined  effect  of  flow  enthalpies  and  heat  generation  in  the 
fuel  cells  determines  the  final  stack  temperature. 

The  effect  of  operating  pressure  on  the  stack  temperature 
is  given  in  Fig.  5.  As  mentioned  above,  increase  in  operating 
pressure  causes  an  increase  in  stack  temperature.  This  makes 
the  fuel  cell  operate  more  efficiently,  which  is  reflected  in 
increased  cell  voltage  and  stack  ac  power  output,  as  shown  in 
Figs.  6  and  7,  respectively. 

The  increase  in  stack  temperature  implicitly  means  that 
there  is  also  an  increase  in  the  exit  flow  temperature  of 
the  stack.  Proportionally,  both  GT  and  PT  would  be  able  to 
produce  more  mechanical  power  through  expansion  pro¬ 
cesses  with  higher  flow  energy.  The  positive  effect  of  flow 
pressure  on  the  PT  ac  power  output  is  demonstrated  in  Fig.  8. 
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Note  that,  the  compressors  and  turbines  are  assumed  to 
operate  at  their  designed  conditions  and  their  respective 
efficiency  is  assigned  with  a  constant  value,  though  they 
should  be  varied  with  operating  conditions.  Provision  has 
been  made  to  use  the  characteristic  maps  to  cater  for  various 
off-design  operating  conditions. 

The  effect  of  flow  pressure  on  the  heat  recovery  is 
presented  in  Fig.  9.  The  heat  recovery  is  referred  to  the 
amount  of  energy  that  is  obtainable  from  recuperator  3, 
which  can  be  used  to  produce  steam  or  hot  water  for 
industrial  and  commercial  applications.  As  expected,  more 
heat  can  be  recovered  from  the  higher  operating  pressure  due 
to  the  higher  flow  temperature.  The  effect  of  flow  pressure  on 
the  system  efficiency  is  shown  in  Fig.  10.  The  simulation 


results  indicate  that  although  the  increase  in  fuel  and  air 
pressures  increases  the  work  required  by  the  compressors, 
the  improvement  in  overall  useful  energy  output  outweighs 
this  extra  work  required,  and  thus  enhances  the  system 
efficiency. 

5.3.  System  performance  under  different  fuel  flow-rates 

The  influence  of  fuel  flow-rate  on  the  stack  temperature, 
cell  voltage,  stack  and  PT  ac  power  outputs  and  system 
efficiency,  is  illustrated  in  Figs.  11-15,  respectively.  In  this 
case  study,  the  pressure  ratios  of  the  fuel  and  air  compressors 
are  fixed  at  9,  the  air  flow-rate  at  I70kmolh  and  fuel 
utilisation  rate  at  80%.  With  a  constant  fuel  utilisation  rate,  a 
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Fig.  15.  Effect  of  fuel  flow-rate  on  system  efficiency. 


higher  rate  of  fuel  flow  means  that  more  chemical  energy  is 
converted  to  electrical  energy.  Accordingly,  more  current 
will  be  produced  and  hence  the  current  density  is  increased. 
This  rise  current/current  density  has  a  linear  relationship 
with  the  amount  of  hydrogen  consumed.  The  effect  of  fuel 
flow-rate  on  stack  temperature  is  presented  in  Fig.  1 1 .  The 
higher  the  amount  of  current  drawn  from  the  SOFC  stack, 
the  higher  will  be  the  overall  overpotential.  This  will 
increase  the  heat  generation  which,  in  turn,  will  raise  the 
stack  temperature. 

The  effect  of  fuel  flow-rate  on  the  cell  voltage  is  given  in 
Fig.  12.  The  combined  effect  of  the  increased  stack  tem¬ 
perature  and  current  drawn  has  a  significant  impact  on  the 
cell  performance.  The  former  causes  a  reduced  overall 


overpotential  while  the  latter  gives  rise  to  an  increased 
overall  overpotential.  In  general,  an  increase  in  fuel  flow- 
rate  at  constant  fuel  utilisation  will  produce  an  initial 
increase  in  cell  voltage  due  to  the  stronger  effect  of  the 
increased  stack  temperature  on  the  overall  overpotential. 
Then,  the  stronger  effect  of  increased  current  density  on 
overall  overpotential  assumes  a  dominant  role  and  causes  the 
cell  voltage  to  fall. 

The  effect  of  fuel  flow-rate  on  stack  ac  power  output  is 
shown  in  Fig.  13.  Intuitively,  an  increase  in  fuel  flow-rate 
causes  an  increase  in  current  and  thus  raises  the  power 
output  of  the  fuel  cell  stack.  The  power  output  due  to  the 
increase  in  current  is  partly  offset  by  the  drop  of  cell  voltage 
at  higher  fuel  flow-rates  as  can  be  seen  in  Fig.  12.  Increase  in 
stack  temperature  and  gas  flow-rate  brings  about  the  advan¬ 
tage  of  increased  turbine  power  output,  and  thereby 
improves  the  turbine  ac  power  output,  as  indicated  in  Fig.  14. 

The  system  efficiency  decreases  with  increasing  the  fuel 
flow-rate,  as  demonstrated  in  Fig.  15.  Though  there  are 
increases  in  stack  and  turbine  ac  power  outputs  as  well  as 
heat  recovery,  the  increase  in  fuel  consumption  causes  a 
decline  in  system  efficiency. 

From  this  study,  it  can  be  seen  that  there  is  always  a  trade¬ 
off  between  the  best  efficiency  and  the  highest  power  output 
of  a  thermodynamic  system.  For  a  specific  SOFC  technol¬ 
ogy,  operation  of  the  stack  at  higher  efficiency  would  mean 
lower  fuel  consumption  or  higher  capital  cost  (a  greater 
number  of  fuel  cells  is  needed).  For  a  system  to  provide 
maximum  power  output,  it  is  necessary  to  determine  the  fuel 
utilisation  rate  which  corresponds  to  the  current  density  that 
will  yield  the  maximum  power  output. 

On  keeping  the  fuel  flow-rate  unchanged,  an  increase  in 
fuel  utilisation  rate  has  a  similar  effect  on  system  perfor¬ 
mance  as  that  of  increasing  the  fuel  flow-rate.  The  increase 
in  fuel  utilisation  rate  would  mean  that  more  fuel  is  being 
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consumed  by  the  SOFC  stack  which  causes  a  higher  current 
density  in  the  cell.  The  cell  voltage  decreases  with  increase 
in  fuel  utilisation,  which  is  attributed  to  the  increase  in 
current  density,  as  can  be  seen  in  a  typical  polarisation  curve 
of  cell  potential  versus  current  density. 

6.  Conclusions 

A  simple,  natural  gas-fed,  hybrid  SOFC  power-generation 
system  has  been  developed.  The  IRSOFC  stack,  which  is  the 
heart  of  the  system,  interacts  with  other  system  components 
and  causes  a  complicated  interdependency  between  them. 
Based  on  this  simulation  study,  the  following  conclusions 
can  be  drawn. 

(i)  With  the  plant  configuration  under  study,  it  is  possible 
for  a  hybrid  SOFC-GT  power  system  to  achieve  an 
electrical  efficiency  greater  than  60%  and  a  system 
efficiency  greater  than  80%. 

(ii)  Given  certain  specific  assumptions,  an  increase  in 
operating  pressure  would  have  some  positive  impact 
on  system  efficiency,  through  improvements  in  cell 
efficiency,  power  output  of  the  PT  and  waste  heat 
recovery. 

(iii)  Keeping  the  fuel  utilisation  unchanged,  the  effect  of 
fuel  flow-rate  on  system  efficiency,  and  vice  versa,  is 
not  favourable.  Though  the  power  outputs  of  the  stack 


and  turbine  show  some  improvement,  the  increase  in 
system  fuel  consumption  outweighs  these  advantages 
and  causes  a  slight  decrease  in  system  efficiency. 
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